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1 Introduction 

A THUNDER actuator is a composite of three thin layers, a metal base, a piezo- 
electric wafer and a metal top cover, bonded together under pressure and at high 
temperature with the LaRC SI polyimid adhesive. When a voltage is applied be- 
tween the metal layers across the PZT the actuator will bend and can generate a 
force. This document develops and describes an analytical model the transduction 
properties of THUNDER actuators. The model development is divided into three 
sections. First, a static model is described that relates internal stresses and strains 
and external displacements to the thermal pre-stress and applied voltage. Second, 
a dynamic energy based model is described that allows calculation of the reso- 
nance frequencies, developed force and electrical input impedance. Finally, a fully 
coupled electro-mechanical transducer model is described. 

The model development proceeds by assuming that both the thermal pre-stress 
and the piezoelectric actuation cause the actuator to deform in a pure bend in a 
single plane. It is useful to think of this as a two step process, the actuator is held 
flat, differential stresses induce a bending moment, the actuator is released and it 
bends. 

The thermal pre-stress is caused by the different amounts that the constituent 
layers shrink due to their different coefficients of thermal expansion. The adhe- 
sive between layers sets at a high temperature and as the actuator cools, the metal 
layers shrink more than the PZT. The PZT layer is put into compression while the 
metal layers are in tension. 

The piezoelectric actuation has a similar effect. An applied voltage causes the 
PZT layer to strain, which in turn strains the two metal layers. If the PZT layer 
expands it will put the metal layers into tension and PZT layer into compression. 

In both cases, if shear force effects are neglected, the actuator assembly will 
experience a uniform in-plane strain. As the materials each have a different elastic 
modulus, different stresses will develop in each layer and these stresses will induce 
a bending moment. When the actuator is released from its flat configuration, the 
differential stresses are relieved as the actuator bends. 


2 Static Model 

The first step in making a model of the THUNDER actuator is to find expressions 
for the relationship between the the force and displacement developed when a DC 
voltage is applied. 

The actuator has three structural layers, a metal base, a PZT middle layer and a 
metal top cover. Fet these three layers be labeled 
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2 . 1 Simple Bending 



Figure 1: The three layers of the actuator with the coordinates used in the model 

c top cover 
b PZT 
a metal base 

A schematic of the layers of the actuator is shown in Figure 1 where the co-ordinate 
system is chosen to be compatible with the IEEE Standard on Piezoelectricity [1]. 
The actuator is initially flat and oriented parallel to the 12 plane with the 3 axis 
perpendicular to the plane of the layers. The 3 axis is also assumed to be the di- 
rection of polarization of the PZT. The bottom surface of the metal base is aligned 
with y = 0. The length of the actuator is /, the width is w, the coordinate of the 
beam in the 1 direction is ~ and in the 3 direction is y. 

Let the thickness of the three layers be h a , In, and In . Introduce the auxiliary 
variables //„ ,///, and y c which are the coordinates of the layer boundaries from the 
lower surface of the metal base, and x ar xb and x c which are the coordinates of the 
centerlines of each layer. 

I/a ha. 

IJb = ha + h b 
lie = ha + hb + h c 
•la. I/a /2 

•l'b = ( IJa T IJb ) /2 

Xc = {lja + IJb + lJc)/‘2 

2.1 Simple Bending 

To fully describe the elastic deformation of each layer we would need to specify the 
three stresses Ti,T 2 and T 3 and (using the indicial notation from the IEEE standard 
[1]) the three shear stresses T 4 , T 5 and T 6 at every location in body of the actuator. 
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2.2 Assumed Strain Model 


However, if we make the assumption that the layers are thin and that they are 
unable to support stresses in the 3 direction, then T 3 = 0. In the simple bending 
model, we also assume that that actuator only bends in the 13 plane and that T 2 = 
0. In reality, both the thermal strain and the piezoelectric actuation will cause both 
bending and stress in the 23 plane. We also neglect the effect of shear stresses, that 
is assume that plane cross-sections of the actuator remain plane under bending, a 
commonly used approximation that holds for moderate shear loads [3]. The one 
dimensional pure bending model gives simple formulas that aid an engineering 
understanding of the mechanisms and parameters of the actuator. An extension 
that includes the out-of-plane stresses is possible if a little involved. 

2.2 Assumed Strain Model 

The strain in all layers due to bending, thermal pre-stress and PZT actuation will 
be a combination of a uniform strain and a bending strain. We will assume that the 
actuator bends with a radius of curvature R or equivalently with a curvature c = 
R~ l and that the deformation is pure bending, i.e. that plane cross-sections remain 
plane. (Positive curvatures are assumed, by convention, to be concave upward and 
so we expect to have negative curvature under the thermal pre-stress.) The beam 
will then have a neutral plane, a distance y 0 from the lowest surface, at which the 
strain is just the uniform strain S 0 . The strain in each layer is then a linear function 
of its distance from the neutral plane. 

S 1 (y) = S 0 -(y-y 0 )/R (1) 

The actuator will bend until the internal stresses are in equilibrium. For pure bend- 
ing, the equations of elasticity are reduced to a single stress 

T\ = YS l 

where Y is the Young's modulus of the material. 

Thermal strains are caused by the different coefficients of thermal expansion of 
the three layers. Let these coefficients be a a , a& and a c . If the temperature at which 
the adhesive layers in the actuator set is 6 degrees above ambient, then the thermal 
strain induced in each layer is 

Si = —ad 

If a potential v is applied between the top metal layer and the metal base, then an 
electrical field is applied across the PZT. In the simple model, this field is assumed 
to be uniform and is given by 

E 3 = -v/h b 
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2.2 Assumed Strain Model 


If the PZT has been polarized in the 3 direction then an electrical field applied in 
the 3 direction induces a strain in the 1 direction given by 

Si = d 31 E 3 = d 3 \V / h b 


where d 3i is a piezoelectric coefficient of the material. 

We now have three sources of strain, bending, thermal and piezoelectric. The 
resulting strains in each layer must combine to give the induced strain of Equation 


h 


.So - 


y - y o 
R 


T a /Y a - ot a 6 0 < y < y a 

< T b /Y b - o b 9 + d 31 E 3 y a < y < y b 
T c /Y c - a c 9 y b <y< y c 


(2) 


These equations can be rearranged to give the stresses induced by the thermal and 
piezoelectric strains for a given radius of curvature R and neutral axis location y 0 , 


T a = Y a (aJ + S 0 - 


y ~ y o 
R 


T b — — d 3 \E 3 + S 0 — 


y ~ y o ■ 

R 


E = Y c (a c 6 + s 0 - 


y - y o , 

R ' 


or if we let y 0 be the location of the plane where ,5'i ( // ) = 0, 

Vo = Uo + RSq 

then 

y ~ vo \ 


Ta = Y a (o a 9 - 


R 


T b — 1 b {ot b 9 — d 3 \E 3 — 


y - yo 

R 


Vi = Y c (a c 9- 


y - y o 
R 


(3) 


(4) 


(5) 


Note that due the uniform strain induced by the thermal and PZT effects, the plane 
of zero strain does not occur at the bending neutral axis y 0 . 

To solve these equations we will need to find the two unknowns, R and y 0 and 
to this end, we will use two equilibrium equations for the net force and net moment 
at a given section of the actuator. 

We consider first the static case, with no inertial forces. At any cross-section of 
the actuator, the forces and moments must be in equilibrium. Thus 


rya ryb rye 

F= T a dy + / T b dy + / T c dy = 0 

Jo Jya J Vb 


(6) 
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2.3 Applied Bending Moment 


rya ryb rye 

M = - T a ycly - / Tyydy - / T c ydy = 0 (7) 

Jo Jya J Vb 

where F and M are the force and bending moment per unit width of the beam. We 
proceed by substituting the expressions for stress and performing these integra- 
tions to get two equations with the unknowns R and y 0 . As the algebra here can 
get rather messy it is wise to proceed in steps, starting with the case of no thermal 
or piezoelectric actuation and an applied moment Mb- 


2.3 Applied Bending Moment 

If a bending moment is Mi, applied at both ends of the beam with 9 = E 3 = 0, (S 0 
is also zero by definition for pure bending) then force equilibrium gives 


rya ryb rye 

-FR = / Y a (y — yo)dy + / Y b {y- y 0 )dy+ Y c {y - y 0 )dy 

■JO Jy a Jyb 


= Y a 

+u 


y 

— - yyo 


+ 1'6 


J 0 


y 

— - yyo 


Vb 


J Va 


r 

2 


yyo 


ye 


= 0 


(8) 


J Vb 


and taking a representative term 


therefore 


and 


y 2 

— - yyo 

yy 22 

Vb ~ Va , ^ 

= 9 (yb y a )yo 



ya 



= h b (x b — yo) 

(9) 

FR = 

1 aF a. J a T ^-b9 b X b \ c h c X c 

— aha + 1'6^6 + l'c^ c )yo = 0 

(10) 

Y a h a x a -j - YhHhXh -j~ Y c h c x c 

(11) 

yo = - 

1 a 9 a + ^bh b + } c fl c 


Expressions like this will be more compact if we introduce the indicial notation, 
where a repeated index i in an expression indicates that the expression is summed 
over all values of i. Thus 


and 


Yhi 


Yh;X; = 


— Y a h a + Yf)h b + Y c h c 

(12) 

Y a h a x a -I - Yj^hj^X}) -I - Y c h c x c 

(13) 
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2.3 Applied Bending Moment 


We can then write an expression for the neutral axis under bending as 


y o 


YihiXi 

) h 


(14) 


Now, if the applied bending moment is Mi, and moments are taken about y = 0, 


M h R = 


rya 


Y a {y - yo)ydy + 


= Y a 


'/ ; //*//< i 


Va 


+ 


(15) 


Again, taking a representative term 


y 3 // 2 //u 


Vb 


Va 


Vb ~ y 3 a Vb - vl 


-yo 


= 4 + h4 - h„.ny o 


(16) 


it can be shown that 


Yh 3 

M b R = -Mi- + ) - f V // .r i// n 


(17) 


The right hand side of this equation is the bending stiffness B of the composite 
beam ( B is the bending stiffness per unit width as Mi, is the bending moment per 
unit width.) Let 

B = 


so that 

More manipulation will clarify the moment equilibrium equation if we introduce 
the auxiliary variables >>'; = — y 0 when. 


+ YihfX? - (YihiXi)y 0 

] l, : ' + v // - - i:) 7/ r :r ’ 


12 


Yhi 


(18) 


B 


YA | 

12 


+ 1 III!' 1 


(19) 
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2.4 Thermal Pre-Stress 


2.4 Thermal Pre-Stress 


Now consider the case of a THUNDER actuator including the thermal pre-stress but 
with no applied voltage. The equations for the three layer stresses reduce to 


T a = Y a {ot a 6 - 
Tb = T&(o&# — 


T c = Y c (a c 9- 


y - y o ■ 
R 

y - yo > 
R ’ 
y - ijo , 
R ’ 


( 20 ) 


We now need to substitute these equations into the force and moment equilibrium 
equations to find R and y 0 The force equilibrium equation is 


o = 


y{a a R0 + yo) ~ — 


+ 


J o 


^ihiOiRO — h 1 ; h ; y ( ) 1 ihiXi 

If this expression is solved for y 0 we have 


a 1 Tr* / 1 / 1 1 / / 1 1 () 1 / h / ( I / 

yo = TT-, = y o — — — R0 


) h 


Yh i 


(21) 


(22) 


To find the static strain induced by the thermal pre-stress, we can apply a bending 
moment to the beam to counteract the thermal bending moment and flatten the 
actuator. The component of the force equilibrium equation due to bending forces 
and the neutral axis then is zero and comparing Equation 22 with Equation 4 we 
have 

V'.A ,AA,. 

(23) 


Y h ot ’ 

Oo,e = TTT u 


Yh, ; 


The negative stress implies that the thermal stress produces a uniform contraction 
of the actuator as expected. 

The moment equilibrium equation is now 


o = i; 


j(a a R0 + y o )-£ 


+ ... 


= ) h x (\ I/O + V// r //n - 

Substitution of y 0 from Equation 22 gives 

Y;h Y;h 

(YihiXiCti - l )R9 = 


12 


— Y;h;X 2 ; 


1 ;//;./ ; . 1 ; h ; O ; ; ( 1 / h ■ ^ ^ ^ 2 (T ;!/;■> ; 


Yh; 


12 


) h 


(24) 


(25) 
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2.5 Piezoelectric Actuation 


or 


B ,,,, Y;h;X;.Y;h;a; , n 

A If, — ' ) ■Ir.r-n — )0 = ) h ir <\ 0 

h i j h ; 


(26) 


Note that the bending stiffness B is independent of the effect of the thermal pre- 
stress. The same expression for Ah, can be found by solving for the bending mo- 
ment required to keep the actuator flat. The moment equation for zero curvature 
is 


M e 


Y n 


'—{ a a Q + S 0 ) 


+ 


J o 


= ) h x a 0 + 17/ r Si (27) 

and when the static thermal strain is substituted we have the same result as before. 


2.5 Piezoelectric Actuation 


The effect of a voltage applied to the PZT layer is to strain that layer which results 
in a compressive stress in the PZT and tensile stresses in the metal layers. By a 
similar reasoning to that used in the previous section, the bending stiffness of the 
actuator will not change under this load. We can calculate the static strain and the 
effective bending moment of the piezoelectric effect directly by assuming that an 
applied moment keeps the actuator flat and by solving for the required moment. 
The stresses in a flat actuator with a piezoelectric actuation are then 


Ta = Y a So 

T b = Y b (S 0 -d 31 E 3 ) 

/ = T Si 

the force equilibrium equation is 

0 = Y.[S 0 y]!; + Y[(S 0 -d 31 E 3 )y]2 + Y[S 0 y^ 

— 1 / h i . s ( i 1 1, h i, 7.3 1 /'S, 

and therefore the uniform piezoelectric strain is 

c , _ A b h b dsiE s 

“ - ~^h~ 

and the piezoelectric bending moment is 


M p 


Y n 


-Sn 


y a 


+ u 


J o 


-(So — d 3 iE 3 ] 


Vb 


+ u 


J Va 


-(So 


Vc 


J Vb 


(28) 


(29) 


( 30 ) 
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2.6 Displacement 




cantilevered 


simply supported 


Figure 2: Cantilevered and Simply Support Actuator Configurations 


= YihiXiSo ~ Y b h b x b d 31 E 3 

= ( ‘ ~ Xb)Ybh b d 31 E 3 

= {yo - x b )Y b h b d 31 E 3 
= -Y' b h b w b d 31 E 3 


(31) 


2.6 Displacement 

There are many possible ways of mounting and loading a THUNDER actuator, each 
of which will result in a slightly different force /displacement characteristic. For 
the purposes of modeling, we assume that the actuator is mounted in one of two 
ways, either simply supported when we will be concerned with displacements at 
the center of the actuator or cantilevered when we will be concerned about dis- 
placements and forces at the free end. Figure 2 illustrates these cases. An advan- 
tage of considering these cases is that the curvature is uniform along the length of 
the actuator under the induced moments of thermal and piezoelectric actuation. 
Other boundary conditions, such as clamped ends, introduce reaction moments 
that make the curvature a function of location. In addition, these reaction mo- 
ments will result in reduced displacements of the actuator and so are undesirable 
in practice. 

The displacement of the actuator is a function of the induced curvature. For 
small displacements we shall assume that this function is linear. THUNDER ac- 
tuators can exhibit large displacements compared to their length and it may be 
necessary to use the nonlinear function for accurate predictions. However, in most 
cases the linear function will give sufficient accuracy. 

The curvature c is given in terms of the actuator displacement u and location ~ 
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2.7 Accounting for the Adhesive Layers 


by 




1 

R 


d 2 u/dz 


(32) 


[1 + (du/ dz) 2 ] 3 / 2 

and for small values the linear approximation is 

d 2 u 
~ dz 2 

The displacement of the cantilevered actuator can then be found as a function of ~ 


(33) 


u(z) = 


2 R 


(34) 


and the displacement at the free end of an actuator of length / is 


UCF = 


l 2 
2R 


For the simply supported actuator 

u(z) = 

and so the displacement at the center is 


Uss = 


z 2 - zl 

2 R 

P_ 

SR 


In order to find the blocked force developed by the actuator, we could apply the 
force necessary to return the actuator to a neutral displacement. Alternatively, we 
can use the principle of virtual work (in the static case) or Hamilton's principle (in 
the dynamic case) to solve for both the force and the electrical parameters, charge 
and current. These models are developed in the next section. 


2.7 Accounting for the Adhesive Layers 

In the development so far we have neglected the effect of the adhesive layers be- 
tween the metal base and cover layers and the PZT layer. In thin THUNDER ac- 
tuators these layers can make up a substantial component of the total thickness. 
Although the adhesive layer is thin and relatively compliant and so does not sup- 
port significant in-plane stress, it does separate the active layers, moving them 
further from the bending neutral axis and so can affect the bending stiffness and 
induced moments. Fortunately, these layers are simple to include in the model. 
Expressions such as Equation 19 for the bending stiffness are evaluated for the five 
layers, metal base, adhesive, PZT adhesive and metal cover rather than just the 
three original layers. This and other expressions can be used assuming summa- 
tion over the appropriate indices for the five layers. 
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3 Dynamic and Energy Based Modeling 

The next step in the model development requires the introduction of dynamics. 
The aim is to find expressions for the resonance frequency, the developed force 
and the electrical input impedance. The most useful approach is to use energy 
based modeling using the principle of virtual work or Hamilton's principle. This 
approach is also more convenient when including effect of more complicated ac- 
tuator geometries as it leads to a quicker formulation of the appropriate relation- 
ships. 


3.1 Energy and the Hamiltonian 


The THUNDER actuator is a coupled electro-mechanical system. As magnetic ef- 
fects can be ignored, Hamilton's principle applied to the actuator can be written 


[ 2 ] 


i-n 


/ [S(T-U + W) + fSu - qS(f>\ dt = 0 (35) 

Jt i 

where U is the mechanical potential energy stored in the actuator, T is the stored 
kinetic energy and W is the stored electrical energy. / is the force applied at a 
location where the displacement is u and q is the charge applied at a location where 
the electrical potential is o. The stored energies are functions of S the full vector of 
strains in a given layer, T the full vector of stresses, E the full vector of electrical 
field strengths, D the full vector of electrical displacements (charge /area) and u 
the vector of local mechanical velocities. 


T 

U 

w 


1 

f puu dV 

(36) 

2 J 
i 

- 1 

V 


f S'TdV 

(37) 

2 J 

V 

1 

2 J 

f E D dV 

V 

(38) 


The prime denotes a transpose and the integrals are taken over the volume of the 
actuator. 


3.2 Stress and Strain 

For a general piezoelectric material the electrical displacement and the mechanical 
stress can be written as linear functions of the applied electrical field and mechan- 
ical strain [1] 

= c E S - e'E 

= eS + e s E 


T 

D 
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3.3 Generalized Variables 


The parameters in this equation are the matrix of material compliances at a fixed 
field c E , the matrix of clamped dielectric constants r s and a matrix of piezoelectric 
constants e. 

For the simple bending model considered so far these equations reduce to func- 
tions of the stress and strain in the 1 direction and the electrical field and displace- 
ment in the 3 direction. 


7 I = US V'/ '.l 

D 3 = t 33 E 3 + 1 SiN 

The first of these equations is familiar from Equation 29. 

3.3 Generalized Variables 

The analysis proceeds by finding expressions for the external variables displace- 
ment, force, charge and voltage and the internal variables N,/i, E 3 and D 3 in terms 
of a set for generalized variables that meet the requirements of the boundary con- 
ditions of a given actuator. If, as is common in practice, we are concerned about 
the low frequency response of the actuator below its first bending resonance fre- 
quency, then a single displacement variable is sufficient. For this case, and with 
cantilevered or simply supported boundaries, a suitable choice of displacement 
variable is the curvature c = /U 1 . The voltage variable is simply the voltage ap- 
plied to the terminals v. We can then write the various internal and external vari- 
ables in terms of c and v. The two independent internal variables are ,5'i and E 3 

Si = c(y — yo) = cy 1 

—v / In, in the PZT layer 
0 elsewhere 

The independent external variables will depend on the choice of boundary condi- 
tions. Let i/( ~ ) be the displacement of the actuator as a function of location. For the 
cantilevered actuator 

u{z) = c = c~ 2 / 2 

and for the simply supported actuator 

u{~) = c $ss(~) = c(z 2 - zl )/ 2 



(43) 

(44) 


(41) 

(42) 


3.4 Transducer Equations 

The transducer equations are a pair of coupled dynamic equations that relate the 
the two electrical variables voltage v and the charge q to the mechanical variables 
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3.4 Transducer Equations 


displacement u and force /. They can be obtained by variational methods by sub- 
stitution of the equations given above into the expressions for energy and taking 
variations by Equation 35. The method, described in more detail in a paper by 
Hagood et al. [2], results in the equations 

Me + Ac — 0 c = $(~ 0 )/ 

0 c + Cv = q (45) 

where M is a mass, I\ is a stiffness, C is a capacitance and 0 is a transduction 
constant each of which is found by integrating energy expressions over the volume 
of the actuator. In contrast to the previous section the following expressions now 
include the width, w, and the area, A = ho, of the actuator. 

I\ is simply a reformulation of the bending stiffness of Equation 19. 

K = A f Yy r2 dy = AB (46) 

Jy 

C is the clamped capacitance 

f s 4 

C = (47) 

lib 

0 is a reformulation of the PZT actuation moment of Equation 31 

0 = / A ^ b f 31 y'dy = A(y 0 - x b )Y b d 31 (48) 

J PZT h 

The mass M and the function 4 ( e (l ) depend on the boundary conditions and the 
location of the applied force. 

M = f p& 2 (z)dV = pihiio f 
Jv Jo 

where p; is the density of the layer i and 4 ( : ) is chosen appropriately for the 
boundary conditions. $ 0 = $(~o) is the value of the appropriate $ at the location 
c 0 at which the force / is applied. For the cantilevered actuator 

Mqf = p;h;iul 5 /20 = nil 4 / 20 

where m is the mass of the actuator, and if the force is applied at the free end 

= l 2 / 2 

For the simply supported actuator 

M ss = pihiwh/VA) = ???/ 4 /120 

and if the force is applied at the center 

$ss(~o) = —l 2 / 8 
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3.5 Resonance Frequencies 


3.5 Resonance Frequencies 

The first resonance frequency of the actuator is given by 

Wi = y/K/M 


and so for the cantilevered actuator 


Wcf 


20 AB 


p'h "7 ' 


4.472 re” 

l 2 V Pihi 


and for the simply supported actuator 


Wss 


120 AB 


i>hirr 


10.95 nr 

l 2 V p t h t 


If the concern is model accuracy around the first resonance frequency of the actu- 
ator, An alternative choice of generalized displacement variable is the amplitude 
of the first mode of vibration. The first mode shape is then be used to calculate the 
integrals for M and K and we obtain the slightly lower first resonance frequencies 
[4] 


wcf 


w S s 


3.52 nr 
l 2 V p t h t 

9.87 nr 
l 2 V Pihi 


(49) 


If the concern is the behavior of the actuator above the first resonance then the 
model should use a set of generalized displacement variables corresponding to 
the amplitudes of the modes of vibration that cover the frequency range desired. 


3.6 Alternate Formulations 

If the excitation is sinusoidal with frequency u and using the relationship u = r <1> 0 
then the transducer equations. Equations 45 can be written 

(A'-w 2 M)«-$ o 0« = $o / (50) 

0 u + 4> ( / ' v = (51) 

In this form of these equations the input variables are / and q and that the output 
variables are u and v. There are four possible alternative forms of the transducer 
equations depending on the choice of inputs and outputs. 
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Inputs Outputs 

force, charge displacement, voltage 

force, voltage displacement, charge 

displacement, charge force, voltage 
displacement, voltage force, charge 

In total there are four possible equations relating each of the three triplets of vari- 
ables ( a. v, /), ( a. v, q), (/, q, a ) and /, q, v). The first two of these are given by Equa- 
tions 50 and 51. Some algebra gives the other two which are 

$lCf + $ 0 Qq = {{K -lo 2 M)C + & 2 ) u (52) 

-$o / + (A' - u 2 M) q = {{K - u 2 M)C + 0 2 ) v (53) 

4 Model Extensions 

4.1 Plate Bending 

The analysis used so far has used a beam model of the actuator. This assumes that 
the effect of the PZT is to bend the actuator in the 13 plane only (Figure 1). In reality, 
both the thermal pre-stress and the PZT actuation cause three dimensional strain, 
with bending in both 13 and 23 planes. A full analysis of this behavior is possible 
but for engineering purposes it is sufficient to modify the beam model by suitable 
approximations to account for these three dimensional effects. 

The first effect of note is perhaps slight. The bending stiffness of equation 19 
was derived assuming the the stress P 2 = 0. Thin plate theory allows for stresses 
in the other plane, results in a bending stiffness modified by the Poisson's ratio a. 
Of course, each material will have a different value of a and the full analysis of the 
multi-layered thin plate will take this into account. However, a good engineering 
approximation is to assume that they are all equal to their average weighed by 
the thickness of the layers. Let a be this average Poisson's ratio, then the bending 
stiffness of the multi-layered plate is 

5 1 plate = (1 _ ( 54 ) 

As a is usually about 0.3, the effect is to raise the stiffness by about 10% and the 
first natural frequency by about 5%. 

4.2 Static Shell Curvature 

Of perhaps greater significance is the effect of static curvature on the bending stiff- 
ness. It is well known that curvature in the 23 plane will increase the stiffness of 
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4.3 Circular Actuators 


a plate in the 13 plane and vice versa. The curvature of the thermal pre-stress has 
the effect of stiffening the actuator for the PZT induced bending, raising the first 
resonance frequency and reducing the displacement. Again, a full analysis of this 
effect is possible, but a useful engineering approximation is fortunately available. 
The theory of shallow spherical shells [6] gives a formula for the increase in natural 
frequency of a curved homogeneous plate as a function of the radius of curvature 
R. 

y 

‘A I shell — u n I plate + ^ 2 

where Y is the young's modulus, p is the density. The equivalent bending stiffness 
is then 

Yhl 4 

^ I shell — B | plate + (^6) 

where the constant [3 2 depends on the choice of boundary conditions. A major 
complication occurs when we wish to calculate the shell bending stiffness due to 
thermal pre-stress as the radius of curvature R itself depends on the bending stiff- 
ness. By equation 26 

R = B/Mg (57) 

where Mg is the thermally induced bending moment. We thus have a cubic equa- 
tion which must be solved for the shell bending stiffness 


B 3 


I shell 


= B 2 


JD I 

shell D I plate 


+ 


Vi, I'M: 

~ 2 


(58) 


4.3 Circular Actuators 

A full analysis of a circular actuator from first principles is possible using a similar 
analysis to that used above for beam and plate actuators. However, as in the pre- 
vious section, it is sufficient for engineering purposes to modify the plate model 
by suitable approximations. 

It was shown in sections 2.4 and 2.5 that the effect of both thermal pre-stress 
and piezoelectric actuation is to induce a bending moment. An equivalent effect 
to the moments induced by distributed action of thermal and piezoelectric strain, 
can be achieved by applying an equal line moment to the edge of the plate. Thus, 
if the rectangular plate is transformed into a circular plate, the line moments at 
the edges are transformed into a line moment around the circumference and the 
formulae of equations 27 and 31 are still valid. As a circular disk certainly does 
not support plane bending, the appropriate bending stiffness is that of equation 
54. The displacement at the center of a simply supported disk of radius a with a 
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line moment M 0 applied on its circumference is [5] 

M 0 ci 2 (l — a 2 ) 
U ° ~ 2B(1 + a) 

and the first resonance frequency is given by [6] 


(59) 


Ui = 



(60) 


The mass term in the transducer model can be found by integrating the square of 
the generalized displacement function 4> 2 {r) over the area of the disk where 


4>{r) = (r — a) (r -\- a) / 2 

and so 

M = f pihi(f>(r)dA = pihia 6 2ir /24 = ma 4 / 12 
J A 

and the other parameters can be readily derived. 


(61) 

(62) 


5 Equivalent Circuit Transducer Model 

It was demonstrated in previous sections that the parameters of the transducer 
model depend on the particular boundary conditions chosen for a given applica- 
tion. In practice, the pure boundary conditions of the analytical model are difficult 
and impractical to implement and the actuator mounting is chosen to meet the me- 
chanical requirements of the application while attempting to optimize one or more 
performance variables. For example, if the THUNDER actuator is to be used as a 
sound source, a mechanical requirement is that a good air seal is obtained from the 
front to the back of the actuator and that this seal does not fail or fatigue during 
operation. Good acoustic performance requires that the volume displacement of 
the actuator at the design frequency be optimized. If the actuator is to be used as 
a source of force, the mounting must allow free movement while transmitting the 
force to the supporting structure. 

Another major practical effect is the static curvature caused by the thermal pre- 
stress which occurs both across the length and width. This shell curvature will 
stiffen the actuator leading to higher resonance frequencies and lower displace- 
ments. The actual curvature of an installed actuator is hard to accurately predict 
as it depends on both the details of the mounting and the manufacturing process 
of the actuator. Small variations in each can lead to large variations in resonance 
frequency. 
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Figure 3: An equivalent circuit for a THUNDER actuator 


Thus while the analytical model is useful for understanding the relationship be- 
tween design variables (material choices and dimensions) and performance vari- 
ables (force and displacement per volt or ampere), it can only provide approximate 
predictions of actual performance variables. An experimentally based model is 
then useful for designing practical actuators and comparing their performance. 

When designing a THUNDER actuator for a particular application, it is useful 
to have a transducer model that relates the electrical input power variables v and 
i to the mechanical input power variables / and ii using parameters that can be 
experimentally measured and hence allow candidate designs to be compared. A 
useful tool for this purpose is the equivalent electrical circuit shown in figure 3. 

In the center of the diagram is an ideal transformer with a 'turns ratio' of N. 
N is the key transduction constant that relates input current to velocity (or charge 
to displacement). Turning the circuit around would give the sensor relationships 
of volts per force and charge per displacement. ( represents the electrical capac- 
itance for the stress free actuator. ( represents the mechanical compliance for a 
open circuit actuator 

Solving for the circuit variables gives the transducer equations in terms of N, 
f t and M m . 


(1 - uj' 2 (C m + N 2 C e )M m ) u - NC e v 
NC e u + C m C e v 
C m f ~\~ N q 

— NC e f + (1 — io 2 (C m + N 2 C e )M m ) q 


in -i- N 2 r \ f 

V w m l 1 y ^ e ) J 

(63) 

{C m + N 2 C e ) q 

(64) 

(1 - id 2 C m M m ) u 

(65) 

C e ( 1 - LO 2 C’m Mm ) V 

(66) 


If the transducer equations derived from the equivalent circuit are compared 
to those derived from the analytical model (Equations 50 to 53) then the following 
equivalences can be found 


N 


0$o 

KC 

.M 

A" 


( l 

\1 + Q 2 /KC 

( l 

VI + & 2 /KC 


(67) 

(68) 
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and 


C e = C{ 1 + & 2 /KC) 
M 


M m = 




The first resonance frequency is given by 




1 


A' 


(C m + N 2 C e )M m M 


(69) 

(70) 


(71) 


which is the resonance frequency for an open circuit actuator q = 0. There is also a 
second slightly higher frequency given by 


2 1 KC + 0 2 

’ 2 “ C m M m ~ MC 


(72) 


which is the resonance frequency for a shorted actuator v = 0. A final but impor- 
tant quantity is the coupling factor k which is a measure of the efficiency of the 
conversion of electrical energy to mechanical energy Assuming that a force and 
voltage are applied such that u and q are both zero. Then by Equation 63 


-NC e v = (C m + N 2 C e )f 

If a small DC charge A q is applied then by Equation 65 

A u = NAq 


(73) 


(74) 


and the ratio of electrical energy supplied to the actuator to the mechanical energy 
supplied by the actuator is 


-/Aii N 2 C e 0 2 

vAq ~ C m + N 2 C e ~ KC + 0 2 


(75) 


6 Design Approximations 

When designing an actuator for a given application, we would like to know the 
effect of particular choices of model parameters on the performance of the actu- 
ator. The equations describing the THUNDER model are moderately complicated 
and the effect of varying a parameter such as the thickness of the base layer on a 
performance variable such as the displacement per volt, is not readily apparent. 

This section discusses a series of approximations or design guidelines that are 
useful when selecting actuator parameters for a particular design. In general, we 
would like to design a transducer that is matched to mechanical impedance of the 
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6. 1 PZT thickness 


medium into which it works. A THUNDER actuator uses the bending process as a 
mechanical amplifier to achieve more displacement than a direct strain PZT device, 
and so displacement performance parameters tend to be of greater importance 
than force performance parameters . 

Assuming then that the objective of the designer is to make a THUNDER device 
with the best displacement output, the most important performance parameters 
are the displacement per coulomb A, the displacement per volt NC e , and the first 
resonant frequency i«i. 

6.1 PZT thickness 

In general, a thinner PZT results in greater actuator displacement at the expense 
of greater actuator force. In a typical application, displacement is of greater im- 
portance that force and so the designer should select the thinnest PZT material 
available. The thickness of available PZT wafers is determined both by the manu- 
facturing process of the wafer itself and by the need to handle the wafer for subse- 
quent processing (as a ceramic PZT is fragile and fractures easily). The thickness of 
the PZT chosen for incorporation into a THUNDER device is thus likely to be fixed 
by its availability. 


6.2 Base Thickness 


The optimum thickness of the metal base layer is a complex function of the selec- 
tion of materials and thickness of the other layers of the actuator. The effect of 
base layer thickness on displacement can be seen by considering its effect on the 
bending stiffness B and the induced piezoelectric moment M p . 


and the curvature is the ratio 


Y'h 3 

B = ‘ 1 + V// ,r J 
12 * 

(76) 

M p = -w b Y' b h b d 31 E 3 

(77) 

c = M p / B 

(78) 


The bending stiffness obviously increases as the base thickness increases. M p also 
increases through the effect of moving the centerline of the PZT .n, further from the 
bending neutral axis y 0 . As a rule of thumb, the maximum curvature occurs when 


Y a h a = Y h h b 


(79) 


and 


'W a = -io b ss h b / 3 


(80) 
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6.3 Actuator Length 


that is the bending neutral axis is just in the PZT (most of the compression strain of 
the PZT is acting on the same side of the neutral axis but no effort is wasted squeez- 
ing base material). With this choice of base thickness, useful approximations of B 
and M p are possible. As PZT has a lower modulus than most metals, the base layer 
thickness will be smaller than the PZT layer thickness and 


B 


Ybtf 

3 


Similarly, 


and so 




— Y h d : 


&M31-E3 


c = 


M p c/3 1 £3 

T* hT 


(81) 


(82) 

(83) 


6.3 Actuator Length 

Having selected the PZT and base layer thicknesses, the next step is to choose the 
length of the actuator. The displacement output increases with the square of the 
length of the actuator and so longer actuators are preferred. However, the first 
resonance frequency of the actuator is inversely proportional to the square of the 
length of the actuator and, in general we will want to operate the actuator at fre- 
quencies below the resonance. At and near the resonance, although the increased 
output is attractive, the performance of the actuator is unpredictable and the phase 
shift accompanying the resonance makes the actuator difficult to control. The de- 
sign choice is then to make the actuator as long as possible such that the operating 
bandwidth of the actuator falls below the first resonance. 

As the thickest layer, the PZT will also dominate the expression for the surface 
density, p t hi ~ p b h b and so the first resonance frequency can be approximated by 


PL^ = 

l 2 V Pihi \f%l 2 V p b 


(84) 


where the constant / 3 depends on the choice of boundary conditions, bearing in 
mind that the stiffness of the actuator will, in practice, depend on the specific de- 
tails of both the mounting and the static curvature induced by the thermal pre- 
stress. 

If loq is the greatest frequency at which the actuator will be operated, then we 
require that 

UJ 0 < LO 1 (85) 
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6.4 Transducer Parameters 


Thus the length of the actuator is given by 

r 2 < 


2 ^ Ph jY b 


\/3iao V p b 


and so the maximum displacement 


,2 , P P m 

«max oc cl oc cl 31 P 3 — W — 

^0 V 66 


(86) 


(87) 


6.4 Transducer Parameters 

Using the approximations for B and iu b given above, approximations for the trans- 
ducer parameters I\ and 0 can be found 


and 


, Y b h\ A 
3 


(88) 

n ^ 33-4 

h b 


(89) 

Y^h^Adsi 

0 sa 
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(90) 

ml 4 

Mcf = 9 q Alss = 

ml 4 

120 

(91) 

9 Yb 

ui,cf — 4 . 6 —— uii'SS — 
1 n b 


(92) 


6.5 Circuit Parameters 

Using the approximations for I\ and 0 given above, approximations for the circuit 
parameters can be found 

N*pp (93) 

h b Ae 33 

The ratio cl 31 / e 33 is the g 31 constant for the PZT material. 

£33 A 
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Cr, 


h b 
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(94) 

(95) 
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6.6 Static Thermal Curvature 


6.6 Static Thermal Curvature 

Accurate estimates of the static thermal curvature are difficult as as the eventual 
curvature depends on both the shell bending stiffness and the precise details of the 
mounting of the actuator. 

An approximation for the thermal bending moment is 

Me = Y b h 2 {a a — a b )9 / 3 

where a a and a b are the coefficients of thermal expansion for the base and PZT 
layers and 9 is the temperature (above ambient) at which the actuator adhesive 
set. An estimate of the curvature is then 


c = Me IB = {a a - a b )9/h b 


If the curvature is large then the actuator will be stiffer than predicted by a 
factor x where 


x 


(l 2 (Q a -Q b )9\ 2/3 

V Hi3 ) 


(97) 


The bending stiffness B and the stiffness I\ must be increased and the compliance 
C m must be decreased by this factor. 


Notation 


B bending stiffness 

C capacitance 

( transduction free capacitance 

C m transduction open circuit compliance 

D electric displacement 

E 3 electric field strength in 3 direction 

F force 

K stiffness 

M mass 

M moment 

Mi, applied bending line moment 

M m transduction mass 

M p piezoelectric actuation line moment 
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6.6 Static Thermal Curvature 


Mfheta ■ ■ 
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. . thermal actuation line moment 

. . transduction constant 

. . radius of curvature 

. . mechanical strain 

. . uniform or static strain 

. . mechanical stress in 1 direction 

. . mechanical stress in 1 direction in layer i 

. . mechanical modulus 

. . mechanical modulus of layer i 

. . radius 

. . curvature 

. . piezoelectric constant 
. . force 

. . thickness of layer i 
. . length 
. . charge 

. . actuator displacement 
. . voltage 
. . width 

. . centerline of layer i 
. . bending neutral axis 
. . upper boundary of layer / 

. . coefficient of thermal expansion of layer i 
. . boundary condition frequency constant 
. . dielectric constant 
. . effective neutral axis 
. . first natural frequency 
. . mode shape 
. . material density of layer i 
. . Poisson's ratio 

. . temperature above ambient at which adhesive set 
. . transduction constant 
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